Strong nonlinear optical response of graphene flakes 
measured by four-wave mixing 



E. Hendry^*, P. J. Hale\ J. J. Moger\ A. K. Savchenko\ and S. A. Mikhailov^ 
1 School of Physics, University of Exeter, EX4 4QL, U.K. 
2 Institute of Physics, University of Augsburg, D-86135 Augsburg, Germany 

We present the first experimental investigation of nonlinear optical properties of graphene flakes. 
We find that at near infrared frequencies a graphene monolayer exhibits a remarkably high third- 
order optical nonlinearity which is practically independent of the wavelengths of incident light. The 
nonlinear optical response can be utilized for imaging purposes, with image contrasts of graphene 
which are orders of magnitude higher than those obtained using Unear microscopy. 



Graphene, a single sheet of carbon atoms in a hexag- 
onal lattice, is the basic building block for all graphitic 
materials. Although it has been known as a theoretical 
concept for some time [1], a layer of graphene has only 
recently been isolated from bulk graphite and deposited 
on a dielectric substrate [5] . The great interest in study- 
ing graphene is driven by its unique band structure and 
many unusual electrical, thermal, and mechanical prop- 
erties [s", ^ . An important stage in the development of 
graphene-based structures was the demonstration that a 
monolayer of atoms can be visible under an optical mi- 
croscope [Hl-lsl. There have also been demonstrations of 
unusual optical properties of graphene: for example, a 
wavelength independent absorption of 2.3% per graphene 
layer [oMllj. Recently, it was predicted that graphene 
should demonstrate strongly nonlinear optical behavior 
at THz frequencies [12], which could lead to important 
device applications. 

Here we present the first experimental studies of non- 
linear optical properties of graphene at visible and near 
infrared frequencies. It is shown that graphene pro- 
duces strong third-order nonlinear optical response, de- 
scribed by nonlinear susceptibility |x^^^| ^ 10~^ esu 
(electrostatic units), which is comparable to that of other 
strongly nonlinear materials, such as carbon nanotubes 
|13l4l7j. In contrast to carbon nanotubes [16], however, 
this nonlinear response is essentially dispersionless over 
the wavelength range in our experiments (emission with 
wavelength of 760 - 840 nm). We also show that the 
nonlinear response of graphene can be exploited for high- 
contrast optical imaging of graphene flakes: the image 
contrasts for monolayer flakes on a dielectric substrate 
are several orders of magnitude higher than those in re- 
flection microscopy [5, 6, 8]. 

To investigate the nonlinear response of graphene, we 
employ the four- wave mixing technique [18]. Figure [TJa) 
illustrates the principle of the method: two incident 
pump laser beams with wavelengths Ai (tunable from 
670 nm to 980 nm) and A2 (1130 nm to 1450 nm) are 
focused collinearly onto a sample and mix together to 
generate a third, coherent beam of wavelength Ag. In 
our experiment, the pump beams are generated by an 
Optical Parametric Oscillator which results in collinear 



6 ps pulses which overlap in time. The incident pump 
pulses are focused onto the sample using a water immer- 
sion objective with a numerical aperture of 1.2, giving 
rise to excitation powers at the sample of ^ 1 mW. 




FIG. 1. (a) Schematic of the experimental layout, indicating 
the pump beams with frequencies cji, 002, and the emission 
beam with frequency uJe- (b) Diagram of energy conservation 
in four- wave process, (c) Band structure of graphene with the 
three photon energies (arrows) involved in four-wave mixing. 

Imaging is achieved by raster scanning of the excita- 
tion beams over the sample and acquiring the emitted 
signal as a function of its position. Images consisting of 
512 by 512 pixels were acquired with a pixel dwell time of 
2.6 /is resulting in an image frame rate of approximately 
0.6 seconds. The nonlinear emission is collected in the 
backward direction using a 750 nm long-pass dichroic 
mirror followed by a 750 nm band pass filter (210 nm 
bandwidth), used to isolate the nonlinear emission from 
the pump beams. The emission is detected by either a 
red-sensitive photomultiplier tube (for imaging) or by a 
spectrometer (see Ref. [19| for more details of the exper- 
imental setup). 

Centro-symmetric materials, such as isolated sheets of 
graphene, do not possess second-order optical nonlinear- 
ities. The third-order nonlinearity of a material is de- 
scribed by the polarization 
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FIG. 2. Emission spectra of a graphene flake excited with 
pump light of different wavelengths, (Ai, A2): (940 nm, 1224 
nm), (950 nm, 1210 nm), (958 nm, 1196 nm), (967 nm, 1183 
nm) and (977 nm, 1168 nm). 

where c^i, 002 and cos are the frequencies of the elec- 
tric field, of the incident beams. In the degener- 
ate four-wave mixing process used in our experiments 
the two mixing frequencies are equal: coi =003. Due 
to energy conservation, the frequency of emitted light is 
cOe = 2a;i — 002^ Figure [H^b). The degree of optical non- 
linearity of a sample is characterized by its third-order 
susceptibility, x*^^^- Dielectric materials (e.g. glasses) ex- 
hibit very low optical nonlinearities, with |x^^^| < 10~^^ 
esu [20]. Conducting materials can have significantly 
higher nonlinear optical responses: for example, free elec- 
trons in narrow-gap III-V semiconductors give rise to 
some of the highest nonlinearities, with |x^^^| ^ 10~^ 
esu |21|. It is also known that large optical nonlinear- 
ities can arise when the mixing or emission frequencies 
(cji, 002 or 00 e) match the energy difference between elec- 
tronic levels [16]. It is interesting to study the magnitude 
and frequency dependence of the nonlinear response of 
graphene, a gapless semiconductor with linear dispersion 
relation. Figure [TJc). 

Graphene and few-layer graphene samples are fabri- 
cated using the standard method of mechanical exfolia- 
tion j2|] and deposited onto a 100 /im thick glass cover slip. 
Prior to investigation in the nonlinear microscope, the 
layer thickness is estimated via contrast measurements 
under an optical microscope using a method similar to 
Ref. [8]. (We calculated the expected reflection coeffi- 
cient of a system with N layers assuming the dielectric 
constant of each layer to be that of graphite.) Figure [2] 
shows the measured nonlinear signal from a monolayer 
flake as a function of emission wavelength Ag for several 
combinations of pump wavelengths Ai and A2. In all 
measurements we observe a clear spike in emission at the 
wavelength corresponding to the condition uJe = 2uji—uj2' 
When the intensity of the pump pulses is varied, the am- 
plitude of the emission peak shows a cubic dependence 
on the intensity, which confirms the third-order nature 



of the response. 

By changing the pump wavelengths, one can investi- 
gate the dispersion of the nonlinear signal. We have 
found that, for equivalent pump intensities, the ampli- 
tude of the emission only changes by ~ 10% over the 
range of incident pump wavelengths. The wavelength 
range that we use is limited at short wavelengths by the 
restriction that coi — 002 must be smaller than the fre- 
quency of the 2D Raman peak in graphene in order to 
avoid coherent ant i- Stokes Raman scattering \2^. On 
the long wavelength side the range is limited by the fil- 
ters used in the experimental set up. We believe that 
the wavelength independence of the nonlinear signal is 
the result of the dispersion relation of graphene, where 
all photon energies in the four-wave mixing process can 
match electronic transitions. Figure [TJc). 

In Figure [3] we compare the images obtained by mea- 
suring the four-wave mixing intensity to those obtained 
on the same flakes in an optical reflection microscope us- 
ing 550 nm light (the standard optical imaging technique 
for thin graphene flakes [1, Q). One can see a striking 
difference in the visibility of the flakes, which are signifi- 
cantly enhanced in the nonlinear image. This difference 
can be quantified in terms of the image contrast, defined 
as 

C=ImL^I^^ (2) 

-^sub 

where Isub and Igr represent the signal arising from the 
substrate and the graphene flake, respectively. As the di- 
electric substrate does not show any measurable nonlin- 
ear signal, this contrast is very large and limited only by 
the noise on the detector: for a monolayer C = 1.7 x 10^ 
compared to C = 0.08 for standard optical reflection mi- 
croscopy [8,]. (Similar enhanced imaging of carbon nan- 
otubes using nonlinear optics has been recently observed 
in Ref. [23']). We expect similarly high contrast for 
graphene on most dielectric substrates, since these mate- 
rials have very weak optical nonlinearity [20] . This makes 
the method of visualizing graphene flakes using nonlinear 
optical response rather universal. 

Using different flakes, we have investigated the behav- 
ior of the nonlinear contrast as the number of graphene 
layers is increased. In Figure [4] we plot the observed 
nonlinear image contrast as a function of the number of 
graphene layers, which is established from reflection mi- 
croscopy [8]. The nonlinear contrast is very sensitive to 
the number of graphene layers, N. It linearly increases 
with increasing the thickness up to A/" ~ 20, while with 
further increase of N the signal decreases towards zero, 
making thick (> 70 A) flakes appear dark in the nonlinear 
images. 

We have modelled this thickness dependence by con- 
sidering a set of A^ parallel graphene layers, each char- 
acterized by the linear-response conductivity cr(a;) and 
the third-order nonlinearity x^^^- Solving the problem 
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FIG. 3. (a) Standard green light reflection images of two 
graphene flakes, (b) Nonlinear optical images measured with 
pump wavelengths of 969 nm and 1179 nm. 



of the transmission of light through such a system we 
get the standard expression for the transmission ampli- 
tude t{u;) = [1 + N27r(j{u;)/c]~'^ and the electric field at 
each graphene layer E = Eotiuo)^ were Eq is the ampli- 
tude of the incident wave. (We assume that the total 
thickness of the A^-layered system is smaller than the 
wavelength of light.) The nonlinear current j2j^-uj2 
duced in each layer by the electric field harmonics coi 
and 002 is then proportional to x*^^^^(Si ^^(^1)^(^2). 
This current generates an electromagnetic wave at the 
frequency 2uJi — U02 with the electric field amplitude 
<^^^_^^^Nt{2uoi — U02) and the intensity 



z:^ -(3) 



^2uJi —002 



Equation (|3]) describes the nonlinear response of the sys- 
tem at the frequency 2cc;i —002- If one assumes that 
the optical conductivity of graphene is dispersionless 
and universal, g{uj) ^ /Ah [1-11|, equation (|3]) gives 
hui-uj2 ^ N'^/{1 + Nire^ /2hcY . This model reproduces 
the trend in the nonlinear thickness dependence on N 
with no fit parameters, Figure [H The saturation and 
subsequent decrease of the nonlinear signal as a function 
of number of layers arises from the decrease of electric 
field inside the material for thicker samples. The differ- 
ence in the exact positions of the maxima in experimental 
and theoretical curves could potentially be accounted for 
by taking into account the finite imaginary part of the 
conductivity. 

To quantify the magnitude of the observed nonlinear 
response in graphene, we compared the signal from mono- 
layer graphene with that of a well characterized, quasi 2D 
material: a thin film of gold evaporated onto a glass mi- 
croscope cover slip. The experiment with the gold film 
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FIG. 4. The contrast in four- wave mixing images as a function 
of the number of graphene layers. 



was performed in the same experimental conditions. It is 
known that a strong nonlinear response of thin gold films 
arises from the excitation of localized plasmons [25|, . 
In Figure in we present the comparison between the non- 
linear emission of a 4 nm thick gold film with that from 
a monolayer graphene flake. One can see that the non- 
linear emission from the gold film is very different from 
that in graphene. For graphene we observe a large peak 
at the emission wavelength, while in gold one can see 
a significantly smaller peak on top of a smooth back- 
ground caused by two-photon luminescence in gold \2^ . 
To compare the magnitudes of third-order susceptibili- 
ties, one has to compare the magnitudes of the peaks, 
Igr and Iau^ with respect to the background. For the 
results presented in Figure [5l the ratio of the peak inten- 
sities in the two materials is approximately ten. From 
this, the nonlinear susceptibility of graphene, Xgr ^ can 



be found from that of the gold film 



10" 



esu 



2 51 . |26[ ) using the relation 



(4) 



where Lgr and Lau is the thickness of the graphene flake 
and the gold film, respectively, and / is the measured 
intensity of the four-wave signal. Given that Lau ~ 10 x 
Lgr, we can estimate that for a single layer of graphene 
esu. 

It is interesting to compare the nonlinear optical prop- 
erties of graphene to those of other carbon materials, 
such as carbon nanotubes 13-17, 28, 29|. Carbon nan- 



otubes can be purely metallic or possess a finite band 
gap, depending on the manner in which the edges of the 
graphene layer join to form the tube. Since nonlinear 
properties are known to be very dependent on the na- 
ture of electronic transitions [l6|, carbon nanotubes ex- 
hibit a broad range of nonlinear polarizabilities in the 
optical and infrared range. For thin films of nanotubes. 
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In particular, four- wave mixing measurements, similar 
to those carried out here, gave 

_ 10-8 esu [17] for 
nanotube films. Due to the singularity in electron den- 
sity of states at the band edges of nanotubes (Van Hove 
singularity) [l^, 0, [29], theoretical predictions suggest 
that homogenous ensembles of identical tubes will ex- 
hibit Ix^^^l which varies over six orders of magnitude on 
varying excitation wavelengths. This strongly resonant 
nature has a typical resonance width AA < 1 nm [l^ . 
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FIG. 5. Emission spectra of a flake excited with pump wave- 
lengths (969 nm, 1179 nm), compared to the emission of a 4 
nm thick gold film under the same experimental conditions. 

While our estimate of the nonlinear susceptibility of 
graphene is comparable in magnitude to that measured 
for nanotubes, we have shown that it is relatively con- 
stant with wavelength. Figure [21 We suggest that this 
frequency independence in the nonlinear response could 
be useful in graphene applications where a constant re- 
sponse is important, for example, in such nonlinear de- 
vices as bistable mirrors for mode-locking, optical lim- 
iters and switchable elements. 

In summary, we have demonstrated that graphene ex- 
hibits a wavelength independent, strongly nonlinear op- 
tical response in the near infrared spectral region, de- 
scribed by a third-order susceptibility |x*^^^| ~ 10~^ esu, 
and shown that nonlinear measurements are a useful tool 
for visualizing graphene fiakes on dielectric substrates. 
Studies of nonlinear properties of graphene represent a 
promising direction for optical experiments on this new 
material. It would be interesting to extend these studies 
to a broad range of frequencies, including the low-energy 
(THz) range. Experimental studies of nonlinear optical 
properties of graphene should be carried out in conjunc- 
tion with theoretical investigations, which are currently 
very limited [12]. 

This work was funded by the RCUK and EP- 
SRC, S.A.M. acknowledges financial support from the 
Deutsche Forschungsgemeinschaft. The authors also 
wish to thank W.L. Barnes, M. Bonn, M.J. Lockyear 
and A.F. Wyatt for helpful discussions. 



[1] P. R. Wallace, Phys. Rev. 71, 622 (1947). 

[2] K. S. Novoselov, et al, Science 306, 666 (2004). 

[3] A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 

(2007); A. K. Geim, Science 324, 1530 (2009). 
[4] A. H. C. Neto, et al, Reviews of Modern Physics 81, 109 

(2009). 

[5] P. Blake, et al., Appl. Phys. Lett. 91, 063124 (2007). 

[6] C. Casiraghi, et a/., Nano Lett. 7, 2711 (2007). 

[7] D. S. L. Abergel, A. Russell, Vladimir 1. Fal'ko, Appl. 

Phys. Lett. 91, 063125 (2007) 
[8] P. E. Gaskell, et a/., Appl. Phys. Lett. 94, 143101 (2009). 
[9] K. F. Mak, et al, Phys. Rev. Lett. 101, 196405 (2008). 
[10] R. R. Nair, P. Blake, A. N. Grigorenko, K. S. Novoselov, 

T. J. Booth, T. Stauber, N. M. R. Peres, and A. K. Geim, 

Science 320, 1308 (2008). 
[11] Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, 

P. Kim, H. L. Stormer, and D. N. Basov, Nature Physics 

4, 532 (2008). 

[12] S. A. Mikhailov, Europhys. Lett. 79, 27002 (2007); S. A. 

Mikhailov and K. Ziegler, J. Phys.: Condens. Matter 20, 

384204 (2008). 
[13] C. Stanciu, et a/., Appl. Phys. Lett. 81, 4064 (2002). 
[14] S. Tatsuura, et al, Adv. Mater. 15, 534 (2003). 
[15] J. S. Lauret, et al, Appl. Phys. Lett. 85, 3572 (2004). 
[16] V. A. Margulis, et a/., Opt. Commun. 249, 339 (2005). 
[17] J. T. Seo, et a/., J. Phys. Conf. Ser. 38, 37 (2006). 
[18] P. Ye and Y. R. Shen, Phys. Rev. A 25, 2183 (1982). 
[19] J. Moger, et al, Opt. Express. 16, 3408 (2008). 
[20] R. Adair, et a/., J. Opt. Soc. Am. B 4, 875 (1987). 
[21] P. A. Wolff and S. Y. Auyang, Semicond. Sci. Technol. 

5, S57 (1990). 

[22] When ooi — C02 corresponds to an optically active phonon 
frequency, one can observe coherent ant i- Stokes Raman 
scattering. Since this paper focuses on the contribution of 
the graphene band structure to the electronic enhance- 
ment of x^^\ excitation of optical phonons is avoided: 
we therefore measure with uji — 002 smaller than the fre- 
quency of the 2D Raman peak (see Phys. Rev. Lett. 97, 
187401 (2006)). 

[23] H. Kim, et al, Nano Letters 9, 2991 (2009). 

[24] D. S. Choi, et a/., Opt. Express 16, 2604 (2008). 

[25] E. Xenogiannopoulou, et aL, Opt. Commun. 275, 217 
(2007). 

[26] D. D. Smith, et al, J. Appl. Phys 86, 6200 (1999). 

[27] H. Kim, et al, Nano Lett. 8, 2373 (2008). 

[28] V. A. Margulis, Physica B 245, 173 (1998). 

[29] V. A. Margulis, J. Phys. Condens. Mat. 11, 3065 (1999). 



